Suppression of Faraday waves in a Bose-Einstein condensate in the presence of an 

optical lattice 
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We study the formation of Faraday waves in an elongated Bose-Einstein condensate in presence of 
a one-dimensional optical lattice, where phonons are parametrically excited by modulating the radial 
confinement of the condensate. For very shallow optical lattices, phonons with a well-defined wave 
vector propagate along the condensate, as in the absence of the lattice, and we observe the formation 
of a Faraday pattern. By increasing the potential depth, the local sound velocity decreases and when 
it equals the condensate local phase velocity, the condensate becomes dynamically unstable and the 
parametric excitation of Faraday waves is suppressed. 

PACS numbers: 03.75.Kk,03.75.Lm 



I. INTRODUCTION 



The term Faraday waves refers to a surface density 
modulation generated by the interference of counterprop- 
agating phonons excited by an oscillatory motion of a 
nonlinear medium [l| . The formation of Faraday waves 
has been studied in several physical contexts Q, includ- 
ing convective fluids, nematic liquid crystals, nonlinear 
optics, biology, and, recently, in Bose-Einstein conden- 
sates (BECs) [3|. In the context of ultracold gases, Fara- 
day patterns can be excited by the modulation of the 
scattering length J4| or by varying the transverse trap 
confinement [3J, |_5|— lFJ , the main ingredient being the mod- 
ulation of the nonlinear interaction. The study of the for- 
mation of spatial structures in a quantum fluid confined 
in a smooth potential, such as the cigar-shaped harmonic 
trap used in the BEC experiment [Sj, can give access 
to information about the excitation spectrum (and thus 
about the equation of state of the system) through the re- 
lation between the excitation frequency and the measured 
wave vector of the Faraday wave. This is somehow similar 
to Bragg spectroscopy on ultracold gases confined in an 
optical lattice, where the excitation spectrum is inferred 
by exposing the system to a periodic lattice modulation 
and measuring the energy absorption [8l-fl2l|. Proceeding 
along this line, one may expect to probe the spectrum of 
an ultracold gas, even in the presence of an optical lattice 
[l3| , by observing the formation (or the lack of formation) 
of a Faraday pattern for any lattice depth provided that 
the system remains compressible (i.e. the lattice does not 
induces a superfluid-Mott transition). 

In this paper we show that the presence of an optical 
lattice may dramatically alter the parametric excitation 
of Faraday waves. In particular, we show that the for- 
mation of the Faraday waves can be suppressed not only 
by the appearance of a gap in the spectrum of collec- 
tive excitations, but also by the fact that the local flow 
of the condensate exceeds the Landau critical velocity 



fl4| . We focus our study on a cigar-shaped condensate 
where phonons are parametrically excited by modulating 
the tight radial confinement at a frequency fl. The para- 
metrically excited phonons have an energy Ml/ 2 and a 
Faraday wavevector qp = q(fl/2) determined by the en- 
ergy spectrum of the system. We add a one-dimensional 
(ID) optical lattice along the axial direction, and, by in- 
creasing the lattice potential depth Vq, we observe the 
suppression of the parametric excitation of phonons with 
wavevector qp, even if Ml/2 is still an allowed energy of 
the collective modes. The sound mode at qp , excited by 
the parametric instability, is submerged by several other 
modes excited by the dynamical instability. This insta- 
bility can be related to the previously-studied instabili- 
ties in bosonic superfluid currents moving with respect 
to an optical lattice potential. In these set-ups the cur- 
rent remains stable if the superfiow momentum does not 
exceed half the recoil momentum [ll| HU. At variance, 
as in our system the condensate center-of mass is at rest, 
the dynamical instability washing out the Faraday pat- 
tern is due to a local supersonic flow through the peri- 
odic potential induced by the continuous driving of the 
radial breathing mode. Indeed, we observe that the crit- 
ical value Vo, C rit at which Faraday waves are suppressed 
corresponds to the critical value of the local BEC flow 



14]. 



The article is organized as follows. In Sec.|TI]we intro- 
duce the nonlinear Schrodinger equations that describes 
the condensate in the elongated geometry. We make use 
of the three-dimensional (3D) Gross-Pitaevskii equation 
(GP) for the evaluation of the low-lying energy spec- 
trum (Sec. IIII[) and of a time-dependent nonpolynomial 
nonlinear Schrodinger equation (NPSE) to numerically 
study the spatial Faraday pattern formation (Sec. IIV|) . 
In Sec. II VI we compare the results for a cylindrical and a 
cigar-shaped trap and show that the suppression of the 
parametric excitations of Faraday waves occurs when the 
local condensate flow reaches the critical value associ- 
ated to the Landau criterion. Our concluding remarks 
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II. THEORETICAL FRAMEWORK 



while the variance is algebraically determined by f(z,t) 
h 



a 2 (z,t) 



■Vi + wnfPF ■ (6) 



We study a 3D Bose condensate trapped into a cigar- 
shaped potential, possibly with a ID periodic potential 
superimposed along the axial direction. The full descrip- 
tion of the system in the limit of zero temperature and 
low density is given by the 3D GP 



2m 



U{v)+gN\^{v,tt 



V>(r,t) 



(1) 

with ^(r, i) the wavefunction of the condensate, m the 
mass of the atoms composing the condensate, N the num- 
ber of particles, and g — 4irh 2 a s /m the interaction cou- 
pling constant with a s the s-wave scattering length be- 
tween particles. In our setup, the trapping potential U (r) 
is the sum of a cigar-shaped harmonic trap plus an opti- 
cal potential V{z) 

U{r) = ^moj 2 ± (x 2 + y 2 ) + ^mu) 2 z z 2 + V(z) , (2) 

with u>± the trapping frequency in the perpendicular 
direction, lo z the trapping frequency in the longitudi- 
nal direction, and uj z 3> wx. The additional poten- 
tial V{z) is either zero or equal to a periodic potential 
V(z) = Vo sin 2 (qBz), with periodicity d = ir/qs- 

Given the above trap geometry, it has been shown in 
Ref. [Hi that a reliable description of the condensate dy- 
namics is given by an effective ID time-dependent NPSE. 
This equation is derived by a variational ansatz for the 
wavefunction 



t/j(r, t) = cj)(x, y, t, a(z, t)) f(z, t) 



(3) 



where the transverse wavefunction </> is modelled by a 
Gaussian function 



(x,y,t,a(z,t)) 
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/na(z, t) 



-(x 2 +y 2 )/2cr(z,t) 2 



(4) 



with a time- and longitudinal-dependent variance a(z,t). 
The validity of this description is based on the assump- 
tion that the transverse wavefunction (f> is slowly vary- 
ing in such a way that the radial velocity can be ne- 
glected. Furthermore, it has been shown that neglecting 
the transverse excitations does not affect the prediction 
for energy and dynamical instabilities [l8|, or the onset 
of Faraday waves for a parametrically excited condensate 
Q. Both the longitudinal wavefunction f(z,t) and the 
variance a(z,t) are determined by the variational princi- 
ple for the energy, and the result is that the longitudinal 
wavefunction is governed by the NPSE 



2m dz 
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muj 2 z z 2 + V(z)- 



l + 3a s N\ff 



^l + 2a s N\f\< 



(5) 



The 3D density profile and velocity field can be obtained 
as 



p(r) = p{z) 



2 / 2 
-T I (7 



v (r) = ^i^MM^W),, (7) 

2mi p 

with p(z) — |/| 2 the integrated ID density. In the fol- 
lowing we will consider the case of a radial confinement 
w_lo/27t = 200 Hz and lattice spacing d = 4680nm; this 
gives a recoil energy E r = ti 2 q 2 3 /2m ~ Q.14hoj± . The 
potential depth Vo will be used as a parameter. 



III. THE LOW-LYING ENERGY SPECTRUM 
WITH PERIODIC POTENTIAL 

We calculate the low-lying energy spectrum of the con- 
densate subject to a cylindrical confinement (u> z = 0) 
with periodic boundary conditions (PBC) for different 
lattice potential depths. Writing Eq. (fT]) in terms of the 
particle density p and velocity v and considering the hy- 
drodynamic limit of the 3D GP we calculate the 
collective modes spectrum u>(q) solving the eigenvalue 
equation 



mu> 2 8p = gV ■ (p V5p), 



(8) 



Po being the equilibrium density in the Thomas-Fermi 
approximation, and Sp(r) = Sp(r±, z)e lqz the density 
perturbation with zero angular momentum propagating 
along the z axis. For that matter we discretize Eq. ([5]) 
in a rectangular domain (r, z), and impose PBC in z, i.e., 
<5p(rj^,z + d) = 5p(r±,z). Figure [1] shows the low-lying 
energy spectrum for the case of No = 3.2 x 10 5 parti- 
cles per site and several values of Vo- As expected, by 
increasing the lattice depth, the lowest-energy band gets 
narrower and an energy gap opens at q = qs- 



IV. THE FARADAY WAVE EXCITATIONS 

In this section we present our numerical results for dy- 
namics of the condensate in the presence of an optical 
lattice. To excite Faraday waves we proceed as follows: 
first we calculate the ground state fo(z) for a static po- 
tential, setting uj±_ = uj±q, then we switch on the trap 
modulation u>x(t) = wj_o(l + e cos Qt) at the frequency O 
and amplitude e. We fix e = 0.1 and choose Q = 2uj±_ hi 
order to be close enough to the natural breathing mode 
in absence of the lattice. 

We focus on the dynamics of the density and veloc- 
ity field of the atoms for increasing values of the lattice 
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FIG. 1: Low-lying energy spectrum in the hydrodynamic ap- 
proximation for the infinity lattice for Vo /E r = (continuous 
line), 50 (dashed line), 130 (dot-dashed line), 190 (circles and 
continuous line), 250 (circles and dashed line), 350 (circles 
and dot-dashed line). 

depth, and compare the case of an infinite cylinder, i.e. 
uj z = 0, and the case of a cigar-shaped BEC, for the 
same number No of particles in the central site. The 
NPSE is numerically solved using a split-step method 
and spatial Fast Fourier transforms (FFT) in a finite do- 
main with PBC. In our numerical calculations, we found 
that is enough to consider a spatial grid with 1024 points 
and a time step 5t ~ 10~ 4 /uj±q to follow the dynamics 
for long times up to several hundreds of the transverse 
frequency period. 

A. The infinite cylinder 

To analyze the formation of a pattern on top of the 
time-modulated density we proceed as follows. At each 
time we separate the part with the spatial periodicity 
of the optical lattice from the total density profile. Ex- 
plicitely, we write p(z, t) = po(z, t) + Sp where 

p (z, t) = A(t) + Bit) cos(2g B z + tp(t)), (9) 

with A, B and tp fitting parameters that characterize the 
renormalized trap, and Sp the remaining fluctuation. The 
Faraday pattern will be analyzed in terms of the be- 
haviour of this fluctuation. 

In Fig. [5] we show the Fourier transform of central 
density p(z = 0) for a lattice depth Vo/E r = 20. At 
ui±flt ~ 300 the oscillation starts departing from a sim- 
ple scaling solution and a different mode gets populated. 
This is seen clearly in the temporal FFT of the central 
density p(z = 0) depicted in Fig. [5] where the frequency 
lu = 0/2 = tu± t o shows up. Shortly after, a well defined 
Faraday pattern is superimposed on the lattice. This is 
shown in Fig. [3] where the linear density p{z) and the 
spatial FFT of its fluctuation 5p(z) are plotted at a time 



FIG. 2: Temporal FFT of the central density from cjx,q£ = 
to 530 with a lattice depth Vo/E r = 20. 



u)±_ ot = 530. It is worthwhile noticing that in the spatial 
FFT of Sp the peaks at q = ±2gg are never present be- 
cause of the choice of the fitting function po in Eq. 
The peaks in the bottom panel at q = ±0.49gs are the 
wavevectors of the counter-propagating phonons giving 
rise to the Faraday pattern, as expected from Fig. [TJ For 
longer times many more modes are populated and a pure 
Faraday excitation cannot be observed. 

For deep lattices (Vo/E r > 230), the frequency £1/2 
lies within the energy gap and no Faraday pattern can 
be formed. However, as Q can intersect an energy band, 
we may expect fluctuations associated to the wavelength 
of the corresponding collective excitation. In particular, 
at Vo/E r = 350 modes at il/2 are not excited anymore, 
while modes at ui = ft are still present. In the spatial 
domain (see Fig. |4]) two main features can be detected: i) 
contributions from higher Brillouin zones at q = ±2qs n, 
stemming from the periodicity of the GS density profile 
which is not contained in the weak Vo limit ansatz for 
po(z,t) in Eq. ([9]); ii) a weak mode in the vicinity of qs, 
originating from the collective mode at uj = (cf. the 
corresponding upper-band in the Fig. [l}. 

Although for intermediate values of Vo, well before ap- 
proaching the edge of the first Brillouin zone, one nat- 
urally expects to excite the pattern, we found that for 
Vo/E r > 50, the mode at uj — il/2 is not seen to increase 
substantially, and therefore it is not possible to identify 
a Faraday pattern. Indeed, we observe several excited 
wavevectors. In order to quantify the spreading in the q- 
space as a function of Vo we calculated the minimum of 
the average of \q\, (|<7|) m m over the excitation Sp during 
a given time interval. This is displayed is Fig. [S] together 
with its dispersion Aq = ((q 2 ) — (l?!) 2 ) 1 ^ 2 (schematized 
as error bars in the same figure) evaluated at the same 
time as (|g|) m i n - For Vo/E r > 50 the dispersion Aq is 
of the same order of qs, more precisely, roughly all q- 
vectors greater than 0.2-0.3 qs are excited. In Fig[5]we 
compare these results with the expected g^-value (solid 
line), as deduced by the low- lying energy spectrum (see 
Fig. [TJ) setting qp = q(fl/2). As already anticipated, 
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FIG. 3: Faraday pattern in a lattice with depth Vo/E r = 20. 
The top panel shows the density profile Sp(z) as a function 
of qBZ at ajj_,oi = 530. The bottom panel shows the spatial 
FFT of the fluctuation 8p(z) at the same time. 
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FIG. 4: The same as in Fig. but for Vo/E r = 350 and 
(jj± t o~t = 150. In this case we observe the contribution from 
higher Brillouin zones 



we can conclude that a clear Faraday pattern cannot be 
identified if Vo/E r > 50, rather before the appearance of 
the energy gap. 



1. The role of the velocity field 

Aiming at understanding the mechanism of suppres- 
sion of the Faraday wave excitation, we compare the lo- 
cal condensate velocity v(z) with the local sound velocity 
c i° c = V ' 9P/ m - At the maxima of the lattice potential, 
the central density p — p/(ir<7 2 ) is minimum, and we ex- 
pect the ratio v(z)/c 1 ° c (z) takes its maximum value. Fig- 
ure |6] shows the behaviour of the maximum value of the 
ratio v(z)/c° c (z) as a function of time, during the para- 
metric excitation, for the case of a shallow (top panel) 
and of a deep (bottom panel) lattice. In the first case, 
where v(z) is always lower than c£ c (z), the increase of 
v{z) corresponds to the onset of the parametric excita- 
tion of Faraday waves, while in the second plot, as soon 
as v(z) becomes greater than c l °°(z) the superfluid be- 
comes unstable according to the Landau criterion, and 
we cannot identify a well-defined Faraday pattern but 
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FIG. 5: Main wavevector (|g|) m in (stars) identified in the den- 
sity fluctuation of the lattice during a time interval Ld± t o~t — 
1500 and its dispersion Aq (schematized as error bars) com- 
pared with the prediction of a Faraday pattern with a wave- 
length qF (solid line), as deduced by the energy spectrum 
evaluated in the hydrodynamic approach. 
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B. Excitations in a confined lattice 

To allow for a more direct application to current ex- 
periments carried out in optical lattices, hereafter we 
consider the effects of the longitudinal harmonic confine- 
ment. For that goal we take ui z /(2ir) = 20 Hz and total 
number of atoms N = 4.7 x 10 6 , parameters that ensure 
the same number of atoms Nq in the central lattice site 
as in the cylindrical confinement. 

In Fig. [7] we show results for Vo/E r = 1. The con- 
tinuous line, that corresponds to po + 55/5, the factor of 
5 having been chosen to zoom the density fluctuation, 
shows the set up of the Faraday pattern at u)±flt = 450 
compared to the ground state at t = (dashed line). 
The presence of the Faraday excitation can be confirmed 
by the inspection of the spatial FFT of the density fluc- 
tuation at z = 0. The lowest \q\ peaks correspond to 
qp = q(Q/2) (see Fig. [1]), while the second main peaks 
correspond to |q| = qp + 2qs, following the lattice peri- 
odicity. The broadening of the peaks with respect to the 
spatial FFT in the cylindrical confinement are due to the 
inhomogeneous density profile. 

Our analysis for the case of a cigar-shaped confine- 
ment is in agreement with our study for the cylindrical 
trap, namely the axial confinement, in the limit where 
u z /lo±o -C 1, does not affect the results shown in the 
previous section. 



FIG. 6: Maximum of the ratio of the local velocity and the 
local sound velocity for an homogeneous BEC attained over 
the lattice as functions of time for Vo/E r = 20 (top panel) 
and 100 (bottom panel). 



several unstable spatial modes. These findings are qual- 
itatively reproduced if one compares the maximum local 
velocity to the sound velocity obtained from the slope of 
the dispersion relationship u)(q). 

Since the velocity of the excitations are proportional 
to the modulation amplitude e, one expects to excite the 
parametric instability without breaking the superfluid if 
e is sufficiently small. However this may provoke a con- 
siderable delay for the onset of the Faraday pattern, in 
both numerical and real experiments. Indeed we have 
investigated this case for e = 0.03 and found that for 
Vq /E r > 100 we do not observe the excitation of the Fara- 
day mode during the modulation up to u)±,ot — 4.5 x 10 4 , 
corresponding to roughly 30s. Since the lifetime of a 
condensate in most experiments does not exceed a few 
seconds, a sizable value of e, as that chosen in our first 
numerical study would be needed to observe the Faraday 
mode. In turn, this would give rise to a fast increase of 
the local flow velocity. 



V. FINAL REMARKS 

The spectrum of a bosonic gas in an optical lattice 
strongly depends on the lattice potential depth. In the 
strongly correlated regime, where there are few parti- 
cles per well, a deep lattice potential induces a quantum 
phase transition towards the incompressible Mott state. 
Otherwise, deep in the superfluid regime where the av- 
erage well occupation is large, by increasing the lattice 
potential depth the Bogoliubov bands get narrower and 
the gaps widen out, as shown in Fig. [TJ The parametric 
excitation of phonons and the observation of the onset of 
Faraday waves is expected to be possible if the excited 
mode is not in the Mott or Bogoliubov gaps 13]. 

In this work we point out another mechanism of sup- 
pression of Faraday waves. During the parametric exci- 
tation, there is a local superfluid flow through the lattice 
barriers. The presence of these barriers do not affect the 
growth of phononic modes with a well defined wave vector 
while the flow velocity remains below the sound velocity 
[Til ]. By increasing the potential depth, this condition 
cannot be fulfilled anymore: the density at the poten- 
tial maxima decreases and thus the local sound velocity 
decreases too. The superfluid becomes unstable and the 
spectrum of the system cannot be inferred by the para- 
metric excitation. 
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FIG. 7: Faraday pattern in an optical lattice with Vo/E r ~ 
1 confined in an elongated trap. The top panel shows the 
"zoomed" density profile po + 56 p at u^^t = 450 (solid line) 
together with the ground-state density profile (dashed line). 
Bottom panel shows the spatial FFT of the fluctuation dp at 
the same time. 
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